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ABSTRACT 

We present results from the analysis of Fe I 630 nm measurements of the quiet Sun taken with the spectropo- 
larimeter of the Hinode satellite. Two data sets with noise levels of 1.2 x 10~^ and 3 x 10""^ are employed. We 
determine the distribution of field strengths and inclinations by inverting the two observations with a Milne- 
Eddington model atmosphere. The inversions show a predominance of weak, highly inclined fields. By means 
of several tests we conclude that these properties cannot be attributed to photon noise effects. To obtain the 
most accurate results, we focus on the 27.4% of the pixels in the second data set that have linear polarization 
amplitudes larger than 4.5 times the noise level. The vector magnetic field derived for these pixels is very pre- 
cise because both circular and linear polarization signals are used simultaneously. The inferred field strength, 
inclination, and filling factor distributions agree with previous results, supporting the idea that internetwork 
fields are weak and very inclined, at least in about one quarter of the area occupied by the internetwork. These 
properties differ from those of network fields. The average magnetic flux density and the mean field strength 
derived from the 27.4% of the field of view with clear linear polarization signals are 16.3 Mx cm~^ and 220 G, 
respectively. The ratio between the average horizontal and vertical components of the field is approximately 
3.1. The internetwork fields do not follow an isotropic distribution of orientations. 

Subject headings: Sun: magnetic fields - Sun: photosphere - Sun: polarimetry 



1. INTRODUCTION 

Only a small fraction of the solar surface is covered by ac- 
tive regions at any one time. The rest is occupied by the quiet 
Sun network and internetwork (IN). Because of their vast ex- 
tension, these regions may contain most of the magnetic flux 
of the solar surface. It is therefore important to determine 
their magnetic properties. 

During the last decade, the analysis of Zeeman- sensitive 
lines has benefited from advances in solar instrumentation and 
spatial resolution. An important contribution has come from 
the analysis of the nearly diffraction-limit ed measurements 
obtained with the spectropolarimeter (SP; 'Lites et al.' '2001') 
of the Solar O ptical Telescope (SOT ; [jsuneta et al. 2008; 
Suematsu efaH 120081: Shimizu et al. 2008|; llchim oto et al. 
2008) aboard Hinode (iKosugi et al. 2 007j) . This instrument 



performs full Stokes polarimetry of the Fe I 630.2 nm lines 
at high angular resolution (0.''32). The Hinode I SV data 
have led to a better determination of the magnetic field 
strength, inclination, and filling factor distributions in the IN 
through Milne-Eddington in versions of the radiative transfer 
equation (Orozco Suarez et al. 2007b|c| ). These resi i lts an d 
the work of "Lites et alj ([2007. .20081). iRezaei et al .' (l2007l), 
Martinez G onzalez et al] (I2008h . iBeck & Rezaeil ([2009), and 



Bommier et al.i (|2009) have helped to settle the main proper- 



ties of quiet Sun IN fields: they are weak (of the order of hG) 
and appear to be very inclined. 

However, there still exists some skepticism about the ac- 
curacy of the magnetic parameters deri ved from ME inver- 
sions of the Hinode m easurements (^e.g. Jde Wijn et al .1120091: 
iBommier et al.l l2009l) . Doubts affect particula r ly the dis- 
tribution of field inclinations. lAsensio RamosI (l2009h and 
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iBorrero & Kobell (1201 Ih . using two different statistical anal- 
yses, recently pointed out that the information contained in 
the Stokes profiles is not enough to constrain the magnetic 
field inclination when the polarization signals are small and 
Stokes Q and U are buried in the noise. As a consequence, 
the true amount of inclined fields may not be recovered accu- 
rately without clear linear polarization signals. The model as- 
sumptions may also induce large uncertainties in the retrieved 
inclinations. 

Always present in real observations, photon noise hides 
the weakest polarization signals and distorts the larger ones. 
There are two ways to partly circumvent this problem: to 
reduce the noise itself and/or to increase the signal. Both 
strategies lead to higher signal-to-noise ratios (S/N), allowing 
for a more precise discrimination between the magnetic field 
strength, field inclination, and filling factor in the quiet pho- 
tosphere. The noise can be reduced by using larger telescopes 
or longer effective exposure times; the signals are enhanced 
by improving the spatial resolution, since this minimizes the 
cancellation of opposite-polarity fields. 

Here we update our results from the inversion of Hinode/SP 
data (Orozco Suarez et al. 2007b, c; Orozco Suarez 2008) by 
analyzing new measurements with significantly better S/N of 
up to 3500. In Section 2 we describe the observations, study 
the amplitudes of the observed polarization signals, and give 
an account of the Milne-Eddington inversion strategy we fol- 
low. In Section 3 the results of the inversion are presented. 
Section 4 discusses the effects of noise on the inferences, tak- 
ing advantage of the availability of two data sets with dif- 
ferent S/N. In Section 5 we determine very precise distri- 
butions of field strength and field inclination for pixels that 
show linear polarization signals well above the noise level. 
The magnetic flux density, the average field strength, and the 
vertical and horizontal components of the field are calculated 
and compared with previous observational determinations and 
magneto-convection simulations in Section 6. To highlight 
the distinct character of IN fields. Section 7 describes the 
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Figure 1. Top: Small IN subfield of 40^^ x40^^ corresponding to data set #1. From left to right: continuum intensity, mean degree of circular polarization, and 
mean degree of linear polarization. The granulation contrast is 7.44%. This area contains a supergranular cell of about 30000 km^. Tickmarks are spaced by 1 
arcsec. Bottom: Quiet Sun region covered by data set #2. The contrast is 7.3%. The x-axis represents time in minutes and the y-axis distance along the slit in 
arcsec. In both cases, the polarization maps have been computed as [T.f=N^ ^(^0//(A/)-E£a^,, y(A/)//(A/)]/(A^r -A^b) and J^^^^ ^yQ(Xi)^ + U(Xi)^/I(Xi)/(Nr: - 
A/b), respectively. N\j is the index where the Stokes V blue lobe of Fe I 630.25 nm starts, and Nr the index where the red lobe ends. Nzc represents the index of 
the Stokes V zero crossing. Nr -A^b = 22 is the number of wavelength points used for the calculations. The circular and linear polarization maps are saturated at 
±0.15% of the quiet Sun continuum intensity. 



magnetic properties of the network as deduced from the in- 
version. Finally, we discuss our results in Section 8. 

2. OBSERVATIONS, NOISE ANALYSIS, AND INVERSION OF 
STOKES PROFILES 

The observations analyzed in this paper come from two 
different data sets taken at disk center. Hereafter they will 
be referred to as normal map observations (set # 1) and 
high S/N time serie s (set #2) iLites et al.l (I2Q07L |2008), 
Orozco Sua rez et al.' ("2007b"!d), 'Asensio Ramos" (12009),^ 



Borrero & Kobel (201 1), and Asensio Ramos & Manso Sain| 



(1201 Ih have also used these measurements. Their main pa- 
rameters are summarized in Table [T] 

In observation #1, the spectrograph slit (of width 0.1 6'0 
was moved across the solar surface in steps of 0.'' 1476 to mea- 
sure the four Stokes profiles of the Fe I 630 nm lines with a 
spectral sampling of 2.15 pm pixel"^ and a exposure time of 
4.8 s. For data set #2, the slit was kept fixed at the same spatial 



location and the Stokes spectra were recorded with an expo- 
sure time of 9.6 s. The completion of data set #1 took about 3 
hours while the time series of data set #2 covers one hour and 
51 minutes. 

The effective exposure time of observation #2 was in- 
creased by averaging seven consecutive 9.6 s measurements. 
This allowed a final exposure time of 67.2 s to be reached, 
which corresponds to a S/N gain by a factor of about 3.8 with 
respect to data set #1. The averaging decreased the spatial 
resolution of the observations to some degree, but the gran- 
ulation pattern is still perfectly visible due to the longer life- 
time of photospheric c onvective structures (about 5 min, see 
iTitle & Schr iivei^'1998') and the excellent stabilization system 
of Hinode (Shimizu et al. 2008). Indeed, the rms continuum 
contrast of the granulation in the high S/N series is only 0.14% 
smaller than in the normal map. This variation is real, since 
the SOT focus position was optimized for the SP observations 
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Table 1 

Log of the observations 



Stokes V 



Data SET #1 
(Normal map) 



Data set #2 
(High S/N time series) 



Date 

Start time (UT) 

FOV 

(Pixels) 
Exposure time 
Stokes V noise (ay) 
Stokes Q,U noise (ctq u) 



2007 March 10 

11:37:37 

302^^x162^^ 

1024x2048 

4.8 s 
1.1 X 10-3 /qs 



1.2 X 10-3 /( 



QS 



2007 February 27 
00:20:00 

1024x727 

67.2 s 

2.9 X 10-4 /qs 

3X10-4/QS 



in the two cases. 

The polarization noise levels (a) are shown in Table[T] They 
were obtained as the mean value of the standard deviation of 
the corresponding Stokes signals in continuum wavelengths. 
Before evaluating the noise, the data were corrected for dark 
current, flat-field, and instrumental crosstalk. The whole pro- 
cess was accomplished using the routine sp_prep . pro in- 
cluded in SolarSoft. The Stokes profiles were normalized to 
the average quiet Sun continuum intensity, /qs, computed us- 
ing all pixels from each data set. 

Figure[T]displays maps of the continuum intensity and mean 
circular and linear polarization degrees for data sets #1 and 
#2. The different nature of the two observations is appar- 
ent. The normal map, as a "snapshot" of the solar surface, 
reveals the spatial distribution of convective cells and mag- 
netic flux concentrations in the quiet Sun. By contrast, data 
set #2 shows the evolution of these structures. For instance, 
granules (intergranules) are seen as bright (dark) horizontal 
streaks. Flux concentrations also produce horizontal streaks 
which last longer in circular polarization. Note that data set 
#2 exhibits many more linear polarization patches than the 
normal map because of the lower noise level. 

Figure [2] shows histograms of the circular (Stokes V) and 
linear (Q, U) polarization amplitudes of Fe I 630.25 nm (top 
and middle panels, respectively) in the two data sets. They 
demonstrate the large occurrence of weak polarization sig- 
nals. The distributions are similar for the high S/N time series 
and the normal map, although some differences exist. For ex- 
ample, the normal map histograms peak at larger amplitudes. 
In both data sets, however, the peaks are close to the corre- 
sponding noise levels: the Stokes V distributions reach their 
maxima at about 1.95crv and 2.35crv, respectively, while the 
linear polarization histograms peak at 2.42crQ,u and 2.33crQ,u- 

To avoid polarization signals that are highly contaminated 
by noise, we only consider pixels with Stokes Q, U (9r V am- 
plitudes larger than 4. 5 cry. This criterion is very restrictive: 
the probability that one pixel showing pure noise in all three 
Stokes parameters be included in the analysis because one of 
the signals exceed by chance the 4. Say level is only 0.19%. 
This number results from the multiplication of three factors: 
the probability of 7 x 10~^ that a single measurement of a zero 
signal exceeds ±4. 5 cry when the noise follows a normal dis- 
tribution, the three Stokes parameters that we consider, and 
the 90 wavelength points used to compute the amplitude of 
the polarization signals. With such a low probability, we can 
be sure that the inverted profiles are real and not due to noise. 

The bottom panel of Figure|2]shows the fraction of the FOV 
with Stokes V signals above a given amplitude. A similar 
curve is provided for the linear polarization signals (Stokes 
Q or U). The vertical lines indicate 4.5 times the noise level 
in Stokes V for data sets #1 and #2. Only 26.0% of the nor- 
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Figure 2. Top: Amplitudes of the observed Stokes V profiles, in units of 
/qs . The thin and thick lines indicate the normal map and the high S/N time 
series, respectively. Middle: Same, for Stokes Q and U (the one with the 
largest amplitude is selected). Bottom: Surface area occupied by signals with 
Stokes V or Stokes Q or U amplitudes larger than a given value. Vertical lines 
mark the 4.5 cr noise thresholds for the two data sets analyzed in this paper. 



mal map area exhibit Stokes V amplitudes above our 4. 5 ay 
noise threshold. This fraction increases to 70.1% in the high 
S/N time series. The area with Stokes Q or U signals larger 
than 4.5crQ,u is smaller, only 2.1% and 27.4%, respectively. 
In summary, the fraction of pixels fulfilling the selection cri- 
terion (i.e., Stokes Q, U, or V amplitudes larger than 4.5 times 
their noise levels) is 26.8% and 72.7% in data sets #1 and 
#2 (corresponding to about 0.7 and 1.6 Mpx). For compar- 
ison, 62.6% and 85.7% of the FOV covered by the two ob- 
servations show Stokes V signals above Say (cf. the 92.6% 
of Martinez Gonzalez et al. 2008 and the 80-90% reported 
by Khomenko et al.M2005,. both from ground-based observa- 
tions). 

In the normal map the selection of IN areas has been done 
manually to avoid the strong flux concentrations of the net- 
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Figure 3. Magnetic field parameters resulting from the inversion of the normal map (top row) and the high S/N time series (bottom row). Shown from left to 
right are magnetic field strengths, inclinations, and azimuths. Pixels in black did not meet the selection criterion and were not inverted. The field strength map is 
saturated at 1 kG so the pixels in white correspond to fields > 1 kG. The FOVs are the same as in Figure 1. Tickmarks are spaced by 1 arcsec in the normal map. 
In the high S/N time series the x-axis represents time in minutes and the y-axis distance along the slit in arcsec. 

as the average Stolces I profile along 1 '^ of the SP slit cen- 
tered on the pixel. The reason is that we cannot perform a 
two-dimensional average since the data is one-dimensional. 
With this approximation we avoid using measurements ac- 
quired more than a minute apart, although the non-magnetic 
profile may not appropriately account for the effects of 
the PSF. Further information about how the PSF cha nges 
the pol arizat ion signals are given in lOro zco Suarez e t al.l 
(2007aFI and ldel Toro Iniesta & Orozco Suar ez (2010). The 
implicati ons of using a local stray ligh t have been ana- 
lyzed by lAsensio Ramos & Manso Sama (1201 1|) . In prac- 
tice, our approach is equivalent to a two-component inver- 
sion in which the non-magnetic atmosphere has only one 
free parameter, (1-/), with / the magnetic filling fac- 



worlc (^see iQrozco Suarez] 120081 for details). In the high S/N 
time series we have just removed the strong magnetic feature 
of positive polarity visible at around l/3rd of the slit (y = 60^0- 
The observations are analyzed using least-squares inver- 
sions based on a one-component, h orizontally homoge- 
neous M E atmosphere, as described by iQrozco Suarez et aP 
(l2007allbll^ . We include a local non-magnetic component 
to account for the reduction of the polarization signals in- 
duced by the instrument's Point Spread Function (PSF). This 
component was called 'loc al" stray-light contamination in 
IQrozco Suarez et aH (12007 aj) . Here we use a different name 
to better reflect the nature of this contribution and distinguish 
it from real stray light, i.e., light that enters the detector be- 
cause of unwanted reflections in the mirror surfaces, sup- 
porting structures, etc. For the normal map, the local non- 
magnetic contribution is evaluated by averaging the Stokes I 
profile within a l''-wide box centered on the pixel. For the 
high S/N time series we take the non-magnetic component 



^ The instrumental PSF considered by these authors included telescope 
diffraction,. CCD pixelation, and SP slit width. Other optical aberrations and 
stray light were not taken into account. 
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Figure 4. Magnetic field strength (left) and inclination (right) probability density functions for IN regions resulting from the inversion of data sets #1 and #2 
(dashed and solid lines, respectively). The PDF for the magnetic field strength corresponding to data set #2 resamples quite well the one corresponding to the 
normal map, despite the very different S/N. These PDFs are derived from the analysis of pixels with Stokes Q, U or V amplitudes larger than 4.5 times the noise 
level. The fraction of the FOV they cover is 26.8% in the normal map and 72.7% in the high S/N time series. Of these pixels, 7.9% and 37.7% show linear 
polarization signals above 4.5crQ^u- 



tor. The two Fe I lines we re inverted simultaneously with 
the MILOS code (Qrozco Suarez & Del Toro Iniestal 120071: 
iQrozco Suarez et al.i i2010ai) . The use of two lines, as 
opposed to only one (e.g., iBommier et al.l E)09), reduces 
the influence o f noise and leads to more accurate results 
(|Qrozco Suarez et al. 2010a ). 

3. INVERSION RESULTS 

Maps of the field strengths, inclinations, and azimuths re- 
trieved from the two data se ts are shown in Figure |3] (se e also 
IQrozco Suarez et al.ll2QQ7bEI and lQrozco SuarezI 120081) . Pix- 
els in black were not inverted because their Stokes Q, U, or 
V amplitudes do not exceed 4. 5 cry. Qualitatively, the inver- 
sions of the two data sets give the same results. Most of the 
fields are weak (of the order of hG), with the stronger field 
concentrations located in intergranular lanes. The fields are 
very inclined, especially in granules. The azimuth is much 
better recovered in the high S/N time series due to the lower 
noise of the linear polarization profiles. 

Figure [4] s hows the probability density function (PDF; 
ISteinerl |2QQ3|) of magnetic field strengths and field inclina- 
tions in the IN for the normal map and the high S/N time 
series (dashed and solid lines, respectively). We recall that 
they were obtained from the analysis of Stokes profiles whose 
amplitudes exceed a polarization threshold of 4. 5 cry. The S/N 
of data set #2 is about 3.8 times better than that of the normal 
map. However, the peaks of the field strength distributions 
are located at nearly the same position (~ 90 G) and have 
very similar widths. A closer inspection shows that the max- 
imum of the field strength distribution for the high S/N time 
series occurs at around 100 G, i.e., at slightly stronger fields. 
These PDFs seem to indicate that the IN consists of hG flux 
concentrations. In general, the field strength PDF does not 
change when using better S/N data, suggesting that the inver- 
sion results are reliable and not affected by the noise of the 
profiles. 

The inclination PDF for the normal map has a maximum at 
90° and decreases toward vertical orientations. Near 0° and 
180° the PDF increases again. The distribution of field incli- 
nations derived from the high S/N data is very similar, with 
some minor differences: first, the amount of nearly horizontal 
fields increases with respect to that in the normal map, indi- 
cating that the smallest polarization signals (not considered 



before) are also associated with highly inclined fields; sec- 
ond, the PDF is narrower, so that the fields are less abundant 
as they become more vertical. Overall, the two distributions 
point to large IN field inclinations. 

4. RELIABILITY OF THE INFERRED MAGNETIC FIELD 
DISTRIBUTIONS 

The polarization signals we measure in the internetwork are 
tiny compared to those found in active regions. As a conse- 
quence, a careful treatment of the noise is important to in- 
terpret them correctly. The first obvious choice to minimize 
the effects of noise was to set a rather conservative threshold 
on the amplitude of the polarization signals, which had to ex- 
ceed 4.5 times the noise level to be included in the analysis. 
We used this criterion to identify and avoid the noisier signals 
that cannot be inverted reliably. 

The threshold of 4. 5 cry was chosen taking into consid era- 
tion the results presen ted by lOrozco Suarez et al.l (12007 ah and 
IQrozco SuarezI (|2008h . There, Hinode/SF normal map obser- 
vations were simulated with the help of 3D magnetoconvec- 
tion models to determine whether it is possible to recover the 
distribution of magnetic field strengths and inclinations with 
a S/N of 1000 and a threshold of 4. 5 ay. The results showed 
that the field strength and the field inclination can be recov- 
ered with mean errors of less than 100 G and 10°, respec- 
tively. For very weak fields (B ^ 100 G), the rms uncertainty 
of the inclination is smaller that 30°, sufficient to distinguish 
between highly inclined and vertical fields. The errors in field 
strength and inclination are expected to be smaller for mea- 
surements with lower noise levels, such as the high S/N time 
series. 

However, th e influ ence of noise is still of concern. 
lAsensio RamosI (2009) analyzed a small region of the nor- 
mal map using Bayesian techniques and Milne-Eddington in- 
versions. His results suggest that while the magnetic field 
strength tends to be well recovered, it is in general not pos- 
sible to constrain the field inclination when the signals are 
very weak and Stokes Q or U do not stand out prominently 
above the noise. A possible consequence of the limited infor- 
mation carried by the Stokes vector when Q and U are buried 
in the noise is an overabundance of horizontal fields. This is 
because inversion al gorithms try to fit the n oise of the linear 
polarization profiles (Khomenko e t al.ll2003l) . 

Also .Borrero & Kobel (201 1) have recently argued that the 
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Figure 5. Magnetic field strength (left) and field inclination (right) distributions for IN regions of data sets #1 and #2 (dashed and solid lines, respectively). They 
have been determined from pixels whose Stokes Q, U or V amplitudes exceed the same absolute threshold of 0.5% the quiet Sun continuum intensity (black 
curves). For the normal map, this amplitude represents 4. Scry (and the corresponding PDFs are identical to those displayed in Figure|4). For the high S/N data, an 
amplitude of 0.5% /qs is 16.4 times larger than the noise. The percentage of pixels fulfilling the selection criterion is 26.8% in the normal map and 16.4% in the 
high S/N data. Of these pixels, 7.9% and 6.8% show linear polarization signals above the noise threshold, respectively. The solid red lines stand for the inversion 
of the same high S/N profiles once their noise has been artificially increased to a level comparable to that found in data set #1. The percentage of analyzed pixels is 
then 16.4%, of which 6.8% show linear polarization signals above 0.5% Iqs. Finally, thegreen lines represent the results of the inversion of data set #1 calculating 
the non-magnetic component as in data set #2. 



noise present in Stokes Q and U can potentially distort the 
Milne-Eddington inferences. These authors showed how in 
some cases inclined fields may be obtained from purely ver- 
tical weak magnetic fields. This happens when the noise in 
Stokes Q and U is interpreted as a real signal. In the weak 
field regime, producing linear polarization in Fe I 630 nm at 
the level of the noise requires large transverse fields. Thus, 
the noise may be compatible with highly inclined fields. An 
additional effe ct is that the recovered fields are stronger than 
the real ones. Borrero & Kobell (|2011) inverted the two data 
sets used in this paper (but only one spectral line, not the two) 
and detected a small shift of the peak of the field strength PDF 
toward weaker fields for the data with better S/N, which they 
interpreted as evidence that noise was playing an important 
role. 

There are reasons to believe that the field strength and field 
inclination distributions we have derived are reliable and not 
significantly affected by the noise, even if most of the pixels 
do not show significant Stokes Q or U signals. The following 
arguments support this conclusion: 

- If noise were affecting the inversions, the fields ob- 
tained from the normal map (the noisier one) would 
be stronger and more inclined than those from the high 
S/N data, as claimed by Borrero & Kobel (2011). Our 
distributions indeed show a small shift when the S/N 
is improved, but in the opposite direction, i.e., toward 
stronger fields (see Figure |4]l. 

- Although only 27.4% of the pixels in the high S/N time 
series have linear polarization signals above 4.5crQ,u, 
the maps of magnetic field inclination and particularly 
the azimuth derived from the inversion are dominated 
by structures which vary smoothly both in time and 
in space (Figure [3]) . These structures are real and not 
caused by the noise. Therefore, the smallest linear po- 
larization signals below 4.5crQ,u also provide informa- 
tion to constrain the field azimuth and inclination. 

- In the absence of clear Stokes Q and U signals it may 
still be possible to gather information about the field 



inclination: first, the noise in Q and U sets limits on 
it, and second, in the weak field regime Stokes I has 
greater sensitivity to the inclination than Stokes Q, U, 
or V (del Toro Iniesta & Orozco Suarez 2010). 

At this point, a solid demonstration that noise is not af- 
fecting the determination of field strengths and inclinations 
would come from the similarity of the results derived from 
pixels which show the same type of polarization signals but 
have very different noise levels. Thus, in Figured we com- 
pare the PDFs calculated from pixels with Stokes Q, U, or V 
signals larger than 0.5% Iqs in both data sets. For the normal 
map (black dashed line), this amplitude threshold represents 
4.5 times the noise level, but for the high S/N observations 
(black solid line) the same absolute threshold corresponds to 
signals 16 times larger than the noise. Unexpectedly, we ob- 
serve some differences in the shapes of the resulting PDFs. 
For example, the peak of the field strength distribution ob- 
tained from the high S/N measurements is slightly shifted to- 
ward stronger fields compared to the normal map. This im- 
plies a smaller abundance of weak fields and a larger fraction 
of strong fields. Also, the two inclination PDFs show a dip at 
around 90°, but less pronounced in the case of the high S/N 
measurements. 

These differences may reflect a better determination of the 
magnetic field parameters from the less noisy signals. Alter- 
natively, they may indicate intrinsic differences between the 
two data sets or problems with the inversion scheme. To iden- 
tify the actual origin of the differences, we have carried out 
two additional tests. In the first one, we have artificially in- 
creased the noise of data set #2 to the level of the normal map. 
Then we have inverted the pixels with Stokes Q, U, or V sig- 
nals above 4.5 times the new noise level. Since this sample of 
profiles has the same noise as data set #1, one should expect 
PDFs similar to those obtained from the normal map (black 
dashed line in Figure O. The differences between the new 
PDFs and the black solid lines should only be ascribed to the 
increased noise. The outcome of this exercise is shown in Fig- 
ure [5] with red solid lines. As can be seen, the new PDFs are 
closer to the ones derived from the high S/N measurements 



Quiet Sun Internetwork Magnetic Fields 




500 1000 
A [mA] 



-500 500 1000 
A [mA] 



Figure 6. Observed (dots) and best-fit (solid) Stokes profiles corresponding 
to a pixel with linear polarization signals above 4.5 times the noise level. 
The field strength, inclination, azimuth, and filling factor retrieved for this 
pixel are 193 G, 104°, 76 °, and 33%. The horizontal lines mark the 4.5a 
threshold. 

(black solid line) than to the normal map PDFs. Compared 
with the original high S/N measurements, the abundance of 
weak fields (0 < 5 < 300 G) decreases and the fraction of 
strong fields (B > 300 G) increases, but the peak of the PDF 
remains at the same position. The fraction of inclined fields 
is also very similar, with slightly less vertical fields. These 
results indicate that noise is not producing the differences be- 
tween the black curves of Figure \5\ since the inversion of 
data set #2 with enhanced noise does not reproduce the nor- 
mal map PDFs. Thus, the cause of the differences must be 
found elsewhere. Another interesting result is the following. 
From the comparison of the PDFs derived from data set #2 
with varying noise levels (bl ack and red solid curves), we do 
not confirm the conclusion of lBorrero & Kobel (201 1) that an 
increase of the noise in Stokes Q and U shifts the PDFs to 
stronger and more horizontal fields. The red curves, much 
more affected by noise than the black curves but otherwise 
coming from exactly the same profiles, do not show any of the 
two effects. This is probably due to the fact that we analyze 
real observations where the fields ar e not completely vertica l, 
contrary to the situation modeled by lBorrero & Kobell (120111) . 
In the second test we examine the influence of the local non- 
magnetic contribution, which is calculated differently in the 
two data sets. To that end, we have repeated the inversion of 
the normal map using a non-magnetic contribution obtained 
as in the case of the high S/N time series, i.e., by averaging 
the Stokes I profiles over V^ along the slit. The results, de- 
picted in Figure [5] with green dashed lines, show that the frac- 
tion of weak fields (below 300 G) decreases and the amount of 
strong fields increases with respect to the dashed black curves. 
The PDF peak shifts to stronger fields. Remarkably, the field 
strength PDF coincides with that calculated from the pixels 
in the high S/N time series that show Stokes Q, U, or V am- 
plitudes larger than 0.5% Iqs. The inclination PDF does not 
change much. Overall, this test suggests that (part of) the dif- 
ferences between the normal map and the high S/N time series 
observed in Figured] may be the result of the different treat- 
ment of the non-magnetic component in the two data sets. 

5. ANALYSIS BASED ON LINEAR POLARIZATION SIGNALS 

To obtain the most accurate results, in the rest of the paper 
we will focus on the analysis of pixels from the high S/N time 



series showing Stokes Q or U signals above 4.5(Jq,u- With 
Q and U clearly above the noise, a very precise determina- 
tion of the field inclination can be made because the inver- 
sion codejisesjh^maxm^ possi- 
ble (As ensio Ram os 2009!: iBorrero & Kobe]|l2011h . The field 
strength is also obtained more accurately, due to the different 
dependence of the linear and circular polarization on B in the 
weak field regime. 

To illustrate how these signals look like. Figure |6] shows 
the Stokes profiles of a typical IN pixel whose linear polariza- 
tion is just above the 4.5crQ,u level. We also include the cor- 
responding best-fit profiles. The fit is quite successful even 
though the signal is close to the threshold and Stokes V is 
slightly asymmetric. Note that the linear polarization sig- 
nals are above the noise threshold (horizontal lines) at several 
wavelength positions, not just one, which allows the inver- 
sion code to determine the field inclination and the azimuth 
accurately because the whole polarization profiles are recog- 
nizable. In this case, the inversion returns a field strength of 
193 G, an inclination of 104°, an azimuth of 76°, and a filling 
factor of 33%. 

An analysis restricted to pixels showing linear polarization 
signals is very reliable because the effects of noise are largely 
suppressed. However, the price to pay is a reduced coverage 
of the internetwork (only 27.4% of the IN surface area) and a 
bias towards inclined and/or strong fields, since those are the 
fields that more easily produce linear polarization signals. We 
do not know exactly the amount of bias, but this possibility 
should be kept in mind. In any case, the data we are consid- 
ering represent more than l/4th of the internetwork, which is 
non negligible. Of all the pixels selected for analysis, 90.3% 
also show Stokes V amplitudes above 4. 5 ay. 

Figure |7] shows our most precise determination of the field 
strength and field inclination distributions in the IN based on 
pixels of the high S/N time series with Stokes Q or U am- 
plitudes exceeding 4.5crQ,u- These PDFs confirm the main 
results derived from the analysis of all the pixels that met the 
inversion threshold, namely that the IN fields are weak and 
very inclined. It is remarkable that the maximum of the field 
strength PDF remains at more or less the same position, with 
only a small shift toward stronger fields (about 130 G). The in- 
clination PDF is narrower than the ones displayed in Figured 
Very likely, this is due to selection effects: the more vertical 
fields are excluded from the analysis because they produce 
smaller Stokes Q and U signals. Therefore, the amount of 
vertical fields can be expected to decrease with respect to that 
found when analyzing pixels with Q, U, or V above 4.5 times 
the noise level. 

Figure [8] shows the filling factors inferred from the inver- 
sion of the profiles in the high S/N time series with Q or U 
signals above 4.5 times the noise level. As can be seen, the 
maximum of the distribution occurs at around 0.25, with a tail 
extending toward filling factors of up to 0.4-0.6. Such val- 
ues are significantly larger than those obtained from ground - 
based observations, and imply that some of the IN fields are 
almost resolved by the Hinode SP. This seems plausible in 
vi ew of the recent d etection of resolved network flux tubes 
by Lagg et al. (2010), using data taken with the IMaX mag- 
netograph (Martinez Pillet et al. 2011) aboard th e SUNRISE 
balloon-borne solar observatory (Sola nki et al.l 2010). The 
resolution of the IMaX measurements is Of' 18, nearly twice 
as good as that provided by the Hinode SP. 

6. AVERAGE MAGNETIC PARAMETERS IN THE QUIET SUN 
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Figure 7. High-precision magnetic field strength and incHnation distributions for IN regions derived from data set #2. The PDFs are based on the inversion of 
pixels whose Stokes Q or U amplitudes exceed 4.5crQu- They represent 27.4% of the surface area covered by the IN. About 90.3% of those pixels show circular 
polarization signals above the noise threshold. 

Our result s partia lly agree with the MHD simulations of 
ISteiner et"an (|2008|) . in which the magnetic field dynamics is 
mainly driven by flu x expulsion and overshooting convection. 
Steiner et al. (2008) computed the horizontally and tempo- 
rally averaged absolute vertical and horizontal magnetic field 
components as functions of height for two different simula- 
tion runs and found a maximum horizontal/vertical field com- 
ponent ratio of about 2.5 at 500 km (see Figure 1 in their pa- 
per). This is comparable with the value we have obtained from 
the inversions. However, the average horizontal and vertical 
magnetic field strengths they find are still smaller than ours. 

Recently, .Danilovic et al.. (.2010.) have presented a set 
of local dynamo sin gulations a n d hav e com pared them 
iLites et al.l 
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Figure 8. Distribution of magnetic filling factors retrieved from the inver- 
sion of the high S/N time series including only pixels whose Stokes Q or U 
amplitudes exceed 4.5crQ u (solid). For comparison, we also represent the 
distribution corresponding to all pixels whose Stokes Q, U, or V amplitudes 
exceed 4. 5 ay (dashed). 

6. 1 . The average field strengths 
We calculated the mean magnetic field strength, {B) = 
Ylf^i^i/^^ ^^^ the mean vertical and horizontal compo- 
nents of the field, (5z) = E!!i I^z,/|/A^=E!!i I^/ cos 7 1 /A^ and 
(^h) = EUBli^Bl,y/yN = Ef.i |5.-sin7|/A^, where A^ is 
the number of pixels, using the results from the analysis of the 
pixels with Stokes Q or U signals above 4.5crQ,u in the high 
S/N time series. The values of (B), {Bj), and {B\^) we obtain 
are 220, 64, and 198 G, respectively. The dominance of (5h) 
over {Bj) indicates that t he fields are highly incli ned in the 
IN, as first pointed out bv iLites et all (I2QQ7L l2QQ8h . Since in 
the rest of pixels the linear polarization signal does not sur- 
pass the noise threshold, they likely have weaker fields. For 
that reason, the above values can be interpreted as upper limits 
for the mean field strength in the internetwork. Note that the 
quantities (5), {Bj) and (5h) are independent of the magnetic 
filling factor. 

The ratio between the horizontal and vertical components 
of the field in IN regions is r ~ 3.1 according to our results. 
iLites et aO (12007'. 2008) estimated a larger ratio r '^ 5 for all 
pixels within the field of view, but this value cannot be directly 
compared to ours because it is based on "apparent" magnetic 
flux densities rather than on int r insic field str engths. In addi- 
tion, the method of ILites et all (I2007L l2008h uses less infor- 
mation than Stokes inversions (for instance, Stokes I was not 
considered) and might be affected by noise differently. 



with the results of iLites et al.l (I2007L l2008i) (see also 
iBello Gonzalez et al.ll2009l for a comparison between ground- 
based observati ons and simulations). Following a similar ap- 
proach as Steiner_eLaD (12008), these authors synthesized the 
Stokes profiles of the 630 nm lines from the simulated mod- 
els and degraded them to the resolution of the Hinode SP 
(iDanilovic et al.ll2008h . After adding noise, they calculated 
the longitudinal and transverse apparent flux densities. Their 
results suggest that current local dynamo simula tions explain 
the value of r -- 5 obtained by Lites et all (l2007L i 2008). How- 
ever, to reproduce the amount of transverse and longitudinal 
flux and the variation of t he flux ratio with heliocentric an- 
gle, IDanilovic et al.l (1201 O f) had to artificially increase the av- 
erage magnetic field strength in the simulation by a factor of 
2 or 3, depending on the dynamo run. With a factor of 3, 
their average field turns out to be abou t 170 G at logr = 0, 
in agr eement with Hanl e measurements (iTrujillo Bueno et all 
'2004"). This prompted I Panilovic et all (1201 0*) to argue that 
current Hinode observations can be well explained by local 
dynamo processes. 

The local dynamo simulations with artificially increased 
fields are roughly compatible with ou r results. The average 
field strength of 170 G worked out by Danilovic et al.l (120101) 
is below the upper limit of 220 G we have deduced. In their 
simulations the vertical and horizontal components of the field 
are also smaller than those reported in this work. Finally, the 
ratio of horizontal to vertical field components varies between 
2 and 4 in the range -2 < logr < -1 (cf. SchUssler & Vogled 
120081) . which is compatible with our value r ^ 3.1 

The two mechanisms put forward to explain the existence 
of very inclined fields in the IN, namely convective overshoot 
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Figure 9. Distribution of magnetic field strength (left) and inclination (right) in network regions resulting from the inversion of data set #1. The solid lines stand 
for all pixels within the selected network areas and the dashed lines for pixels associated with network flux tubes (see text). 



(ISteiner et al .1120081) and local dynamo action toanilovic et al.l 
I201QI) . may operate simultaneously or may not occur in the 
real Sun. To distinguish between the different possibilities it 
is necessary to perform detailed comparisons of the magnetic 
field distributions predicted by the simulations and those ob- 
tained from the observations. This would give much more 
information than just a single parameter such as the ratio of 
horizontal to vertical magnetic field components. 

6.2. The average magnetic flux density 

The determination of the average flux density of IN regions 
has been pursued for many years, which has resulted in more 
than forty papers to date. Unfortunately, there is a large dis- 
parity between the values obtained by the different authors. 
One of the reason is that the estimates are strongly affected by 
the angular resolution of the observations. Also the different 
analysis techniques have contributed to the discrepancies. As 
a result, the flux valu es reported in the literature v ary from the 
2-3 Mx cm-2 of, e.g. JLin & Rimmeld (119991) andlKelleret al.l 
(119941) , to the 1 5-20 Mx cm~^ found by Dominguez CerdefiaJ 



2003 ). IViticchie& Sanchez Almeida (2011), B eck & Rezaeil 
( 20091) . and |BeUo Gonzalez & Kneer (2008). Studies us- 
ing simultaneous observati ons of visible and infrared lines 
give 11-15 Mx cm"^ (Kh omenko et al.l [2005h . An upper 
limit to the flu x density seems to be '^ 50 Mx cm~^ (e.g., 
lFauroberterani2QQl.) . 

We have calculated the average magnetic flux density ((/)) = 

X]/=i 1 4'i I /^ using the parameters retrieved from the inversion 
of the high S/N time series. The flux of individual pixels is 
computed as (/) = fB cosj. Note that, as opposed to the mag- 
netic field strength B, the flux density depends on the mag- 
netic filling factor. The average is carried out in two different 
ways: taking into account only those pixels with Stokes Q or 
U amplitudes above 4.5 ctq^u and considering all pixels. In the 
latter case, the average flux represents a lower limit because 
we assign zero fluxes to pixels without clear linear polariza- 
tion signals (even if they show large Stokes V amplitudes). 

We find a flux density of 16.3 Mx cm"^ in the 27.4% of the 
FOV with significant linear polarization signals. A lower limit 
for the flux, considering all pixels, is 4.5 Mx cm~^. These 
res ults agree w ith the latest estimations of 11-16 Mx cm~^ 
by iLitesI (1201 lb . Our flux densities are slight smaller than 
those derived recently from the analysis of near infrared lines 
by Beck&Rezaei (:2009j), who found fluxes of about 26 G 



and a lower limit of 22 G. They are also smaller than the 
ones obtained from observations in the visible with the GRE- 
GOR Fabry-Perot spectrometer attached to the German Vac- 
uum T ower Telescope: according to Bello Gonzalez & Kneed 
(I2008h . a lower limit for the flux would be 17 G. Finally, the 
net flux density, calculated as (^)n = Xl/^i (l^i/N, amounts to 
-2.9Mxcm"^. 

7. DIFFERENCES WITH THE NETWORK 

In this Section we show that the properties of network and 
internetwork fields are rather different, as expected from their 
different nature. The distribution of network field strengths 
and inclinations obtained from the analysis of the normal map 
can be seen in Figure [9l The PDF for the field strength (solid 
line) peaks at '^ 100 G but has a prominent hump centered 
at about 1 .4 kG, in contrast to what is found in the IN where 
the majority of fields are in the hG range (Figure |7). In the 
same figure we give the distribution of fields for the inner pix- 
els of the network flux concentrations (dashed line). They 
are selected by their circular polarization amplitudes, which 
have to be larger than 0.03 Iqs. For those pixels, the peak at 
100 G shifts toward stronger fields and becomes smaller than 
the hump at 1.4 kG. The PDF vanishes very rapidly below 
300 G. 

The inclination PDF corresponding to the network areas 
(solid line) is dominated by inclined fields, but has secondary 
peaks at ~ 10° and 170°. The slightly asymmetric shape of 
the distribution is due to the predominance of positive polar- 
ities among the network areas selected for analysis. When 
we consider only the stronger network concentrations (dashed 
line), the peaks at 10° and 170° become dominant and the in- 
clined fields disappear. 

The network PDFs show what is expected from intense 
field concentrations, i.e., kG fields with predominantly ver- 
tical orientations. The field strength and field inclination 
peaks we observe are compatible w ith the results of an aly- 
ses at lower spatial resolution (^e.^.. [ Solanki et"an Il987l and 
ISanchez Almeida & Martinez Pilletlll99 4'). The large fraction 
of inclined fields as well as the tail and hump of the field 
strength distribution toward weak fields represent the contri- 
bution o f network canopies t o the PDFs. Network canopies 
(see e.g. iPietarila et al .1120101 and references therein) are as- 
sociated with more inclined and weaker fields and can be 
observed at the vicinity of strong network flux concentra- 
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Figure 10. Small area of about 5^^ x5^^ showing the degree of circular polarization (left), the field strength (middle) and the field inclination (right) resulting 
from the inversion of normal map data. This area contains several network elements. The white pixels were not inverted. Contour lines enclose areas with field 
strengths greater than 700 G (dotted), inclinations lower than 25° (solid), and polarization signal amplitudes above 2 x 10~^/qs (dashed). 



tions. To illustrate this, Figure [TO] displays the magnetic field 
strength and inclination inferred from the ME inversion of a 
small area in data set #1 that includes network elements. Note 
that pixels corresponding to the innermost part of the network 
elements show strong vertical fields, while those at the pe- 
riphery (canopy areas) are associated with weak fields and 
moderate-to-large inclinations. 

8. DISCUSSION AND CONCLUSIONS 

In this paper we have carried out a Milne-Eddington anal- 
ysis of high spatial resolution measurements of the quiet Sun 
taken with the Hinode spectropolarimeter. To infer the mag- 
netic field vector we have applied the i nversion strategy de- 
scribed by lOrozco Suarez et aH (l2QQ7ab to two observations 
with different noise levels: 1.2 x 10"-^ and 3 x 10~'^/qs. 

The magnetic field distributions obtained from the two data 
sets are rather similar, suggesting that they are not biased by 
photon noise. A more detailed analysis reveals that differ- 
ences exist between the distributions when we consider the 
same type of polarization signals (Stokes Q, U, or V larger 
than 0.5% Iqs) but different noise levels. We performed two 
tests to assess how the noise and the calculation of the non- 
magnetic component affect the inversion results. We reach the 
conclusion that noise is of little concern. Rather, the results 
may be influenced by the different method we use to obtain 
the non-magnetic component in the two data sets. This test 
suggest that the field strengths may be slightly over-estimated 
in the high S/N data. The field inclinations do not seem to 
depend on the exact way the non-magnetic component is cal- 
culated. 

We also provide magnetic field distributions in the IN based 
only on pixels that show significant linear polarization signals 
(Figure |7]). Most of them also have large circular polarization 
signals. By restricting the analysis to those pixels we make 
sure that the inversion results are accurate. The downside is 
that they represent only 27.4% of the surface area covered 
by the IN. Our choice also implies that these high-precision 
magnetic field distributions may be biased toward stronger 
and/or more inclined fields — the ones that produce the larger 
linear polarization signals. A better coverage of the IN re- 
quires more sensitive measurements that are not yet available. 

Keeping in mind these considerations, we discuss the main 
results of our analysis in the following subsections. 

8.1. Internetwork fields are weak 



The inversion of Hinode/SF data presented here indicates 
that IN fields are weak, at least in the 27.4% of the FOV 
showing Stokes Q or U signals well above the noise. This 
is in qualitative agreement with the picture derived f rom 
the m ore magnetically sensiti ve Fe I lines at 1565 nm (JLinl 
[i^;lLin & Rimmele 1999; CoUadosll2QQll:lKhome nko etaD 
|2003; Dommguez Cerdena et al. 2006; Beck & Rezaeil l2009h 
and with the simultaneous inversion of the Fe I 630 nm and 
1565 n m lines performed by [Martinez Gonzalez et al. (20^. 
Also Rezaeietail (l2007h and Beck & Rezaei (2009) found 
weak fields in the IN from the analysis of Fe I 630 nm ob- 
servations take n with the Polarim etric Littrow Spectrograph 
( Schmi dt et^ 120031 : iBeck et al.l 12005.) at the German Vac- 
uum Tower Telescope in Teneri fe. Our findings are compati- 
ble with the results obtained by lAsensio RamosI (|2009|) from 
the same data using Milne-Eddington inversions. 

There are other measurements that point to the existence 
of weak fiel ds in the IN. For ex ample, the Mn I 553 nm line 
analyzed by iLopez Ariste et al.l (|2006) suggests that the IN 
is dominate d by fields below 600 G. Us ing the Mn I line at 
1526.2 nm, lAsensio Ramos et aP (l2007h found a Gaussian- 
shaped distribution of field strengths centered at around 250- 
350 G. The manganese lines are important because the hyper- 
fine effects they show make it possible to derive the strength 
of the field directly from the shape of the profiles. Since 
the polarization amplitudes are not used, the inferred field 
strengths are free from errors due to uncertainties in the ac- 
tual stray-light contamination or in the magnetic filling factor. 

The field strength PDF displayed in Figure |7] is not mono- 
tonic and has a clear maximum in the region of weak fields. 
Although the exact location of this maximum is still un- 
der debate, the peak itself is most probably solar in origin. 
Several authors have argued that peaks may result from the 
cancellation of magnetic flux at sub-resolution scales (e.g., 
Martinez Gonzalez et al. 2008), but we do not favor this in- 
terpretation because our spatial resolution is much better than 
that of any previous measurement, which decreases the possi- 
bility of cancellations. This does not mean that cancellations 
do not exist at 0.''3, but that they are less frequent. 

An important question is whether we should have con- 
sidered more complex magnetic topologies than those al- 
lowed by the Milne-Eddington approximation. I Viticchie et al.l 
(1201 1|) have recently analyzed Hinode/SF data assuming 
that the magnetic fields of the IN are small fibrils with 
sizes below ~ 100 km, the mean photon free path in 
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Figure 11. Magnetic field inclination PDF of IN regions in the high S/N 
time series where the linear or circular polarization signals are larger than 
1.35 X 10~^ Iqs (corresponding to 4. Scry, solid line) and pixels where the 
linear polarization exceeds the same threshold (dashed line). The dotted line 
represents the distribution of inclinations for magnetic field vectors with ran- 
dom orientation (isotropic case). Its shape is given by a sine function, here 
represented in logarithmic scale. 

the solar photosphere. More specifically, they used a 
Micro- Structured Mag netic Atmosphere (MISMA ) of the 
type describ ed by Sanchez Almeida et ajj (119961) . With 
this model, IViticc hie et al. (2011) were able to fit the 
asymmetries of the Stokes p rofiles recorded by Hinode 
(|Viticchie & Sanchez Almeidal 12011^) . They found that a 
broad range of field strengths (from hG to kG) are present in 
IN regions. In particular, kG fields would dominate deep pho- 
tospheric layers and hG fields the layers above. They argued 
that an IN consisting of fibrils with kG fields is partly com- 
patible with the field strength PDFs shown in Figure lU be- 
cause the ME inferences carry information from high photo- 
spheric layers. However, we doubt that our findings and their 
results are truly compatible. First, the ME inferences cannot 
be assigned to specific atmospheric layers but rather provide 
average values over the region of formation of the spectral 
lines. Therefore, our field strength distributions should be 
similar to the ones they obtain, w hich is not the case. Sec- 
ond, the inversions performed by IViticchie et al.l (120111) are 
very dependent on the assumptions of the MISMA model. In 
the absence of compelling proofs that the quiet Sun IN con- 
sists of small-scale, mixed-polarity magnetic fibrils that are 
still not resolved by Hinode, we tend to prefer more simplis- 
tic approaches like the Milne-Eddington inversions described 
in this paper. 

Finally, we po i nt out that our results seem to agree with 
those of Stenflo (2010). He reported the existence of two 
kinds of magnetic structures in the quiet Sun, one character- 
ized by strong and vertical fields, and another by weak and 
isotropic fields. We believe that the strong fields are associ- 
ated with the network ele ments, while the weak ones are those 
found in the IN, although lStenflol ([2010) cautioned that the as- 
sociation of IN fields with the weak fields is not completely 
clear. 

8.2. Internetwork fields are highly inclined 

The distribution of field inclinations deduced from the 
27.4% of the FOV with clear Stokes Q or U signals in the high 
S/N time series suggests that IN fields are very inclined. This 
result may be biased by the requirement of significant linear 
polarization, but it is at least valid for one quarter of the IN. 
The existence of nearly horizontal hG fields in IN regions is 
compatible with the large transverse magnetic fluxes found by 



iLites et al.l (I2007L l2008h and lMartinez Gonzalez et aH (l2008h . 
We caution, however, that "very inclined" does not mean 
"horizontal": horizontal fields imply inclinations of 90°, 
which are indeed the most abundant but not the only ones. 

An aspect that is being discussed intensely in the liter- 
ature is whether or not the inclination PDF is compatible 
with an isotropic distribution of magnetic field vectors in 
the FO V dAsensio Ramos"2009'; "Bommier et al."2009|; iStenflol 
I2OIOI: see [Sdnchez Almeida & Martinez Gonzalez 201JJ for a 
review). Based on the shape of the PDF, several authors have 
in fact suggested that IN fields are "turbulent". Figure [TT] 
shows the distribution of magnetic field inclinations deduced 
from pixels in the high S/N time series with circular or lin- 
ear polarization signals above 4.5 times the noise level, corre- 
sponding to amplitudes of at least 1.35 x 10~^ Iqs (solid line). 
We also represent the PDFs of a random (i.e., isotropic) dis- 
tribution of magnetic fields (dotted line) and of pixels in the 
high S/N time series that show linear polarization amplitudes 
above 1.35 x 10~^ Iqs (dashed line). As can be seen, neither 
of the two distributions derived from the Hinode data are com- 
patible with an isotropic or quasi-isotropic distribution of field 
vectors. The amount of very inclined fields obtained from the 
inversions clearly exceeds that of a random distribution. Thus, 
we conclude that IN fields are probably not isotropic, at least 
those that show prominent linear polarization signals (repre- 
senting 27.4% of the solar IN). A detailed analysis by Borrero 
& Kobel (201 1) also points in the same direction. 
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